UPPER BOUNDS FOR THE DENSITY OF SOLUTIONS OF 
STOCHASTIC DIFFERENTIAL EQUATIONS DRIVEN BY 
FRACTIONAL BROWNIAN MOTIONS 



FABRICE BAUDOIN, CHENG OUYANG, AND SAMY TINDEL 



Abstract. In this paper we study upper bounds for the density of solution of stochastic 
differential equations driven by a fractional Brownian motion with Hurst parameter 
H > 1/3. We show that under some geometric conditions, in the regular case H > 1/2, 
the density of the solution satisfy the log-Sobolev inequality, the Gaussian concentration 
inequality and admits an upper Gaussian bound. In the rough case iJ > 1/3 and under 
the same geometric conditions, we show that the density of the solution is smooth and 
admits an upper sub-Gaussian bound. 
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1. Introduction 

Let B = {B^,...,B'^) be a d dimensional fractional Brownian motion (fBm in the 
sequel) defined on a complete probability space {Q,J^,¥), with Hurst parameter H G 
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(0, 1). Recall that it means that i? is a centered Gaussian process indexed by IR+, whose 
coordinates are independent and satisfy 

(1) E[(s/-5^y] = for s,teR+. 

In particular, by considering the family {B^; H G (0,1)}, one obtains some Gaussian 
processes with any prescribed Holder regularity, while fulfilling some intuitive scaling 
properties. This converts fBm into the most natural generalization of Brownian motion 
to this day. 

We are concerned here with the following class of equations driven by B\ 

(2) Xt = x+ f Vo{X:)ds + Y,f V.iX:)dBl, 

Jo ^0 

where x is a generic initial condition and {Vi] < i < d} is a collection of smooth vector 
fields of M*^. Owing to the fact that fBm is a natural generalization of Brownian motion, 
this kind of model is often used by practitioners in different contexts, among which we 
would like to highlight recent sophisticated models in Biophysics [201 ISH l30] . 

As far as mathematical results are concerned, equation ([2]) is now a fairly well under- 
stood object: existence and uniqueness results are obtained for if > ^ thanks to Young 
integral type tools [32||27]. while rough paths methods [T3||22] are required for j < H < ^ 



4 ^ ^ 2" 

Numerical schemes can be implemented for this kind of systems [TTl [13], and a notion of 
ergodicity is also available [T6 | fT7 ] . Finally, the law of has been analyzed by means of 
semi-group type methods [UlSl] and its density has also been investigated in [2| [71 [T9| |28] . 

In spite of these advances, concentrations results and Gaussian bounds for the solution 
to ([2]) are scarce: we are only aware of the large deviation results [23] in this line of 
investigation. The current article is thus an attempt to make a step in this direction, by 
analyzing a special but nontrivial situation. 

Indeed, we consider here equation driven by a fBm with Hurst parameter H G 
(|, 1), and we suppose that our vector fields Vq, . . . , fulfill either of the following non- 
degeneracy and antisymmetric hypothesis: 

Hypothesis 1.1. The vector fields Vq, . . . ,Vd are -bounded, and Vi, . . . ,Vd satisfy 

(i) For every x G W^, the vectors ■ ■ ■ , Vd{x) form a basis ofW^. 

(ii) There exist smooth and bounded functions ufj such that: 



(3) [V,,V,]=J24^k, and ^ = -oo: 



_ J 

ij ~ "^ik- 

k=l 



The second assumption (ii) is of geometric nature and actually means that the Levi- 
Civita connection associated with the Riemannian structure given by the vector fields V^'s 
is 

VxY = \[X,Y]. 

In a Lie group structure, this is equivalent to the fact that the Lie algebra is of compact 
type, or in other words that the adjoint representation is unitary. Such geometric assump- 
tion already appeared in the work [3] where it was used to prove a small-time asymptotics 
of the density. 
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Hypothesis 1.2. The Hypothesis is satisfied and moreover, the vector fields Vi, Vd 
form moreover a uniform elliptic system. That is 

\v'^VV^v\ > A|i;p, for all v G W^. 

Here V = (V^*)ij=i,.--,rf X is a positive constant. 

When H > ^, under Hypothesis 11.11 our main result can be loosely summarized as 
follows (see Theorem 13 . 1 41 for a precise statement): 

Theorem 1.3. Fix H > ^. Let be the solution to equation and suppose As- 
sumption \l.l\ is satisfied. Then for any t G M.*^, the random variable Xf admits a smooth 

density px{t,-). Furthermore, there exist 3 positive constants c[^\ c[^\ c^'^l such that 

Vxit.y) < c^^^ exp ^-cf^ {\y\ - eg) ^ , 

for any y . 

We don't claim any optimality in the quantities c^^\cf''^ and cfl above (whose exact 
definitions are postponed to Section IX^ . Nevertheless, this is (to the best of our knowl- 
edge) the first Gaussian type bound available for solutions of differential equations driven 
by fBm. 

Let us say a few words about the strategy we have followed in order to prove Theo- 
rem 11.31 It is mostly based on stochastic analysis tools, and particularly on a general 
integration by parts formula giving an exact expression for the density px{t^ ■) in terms 
of Malliavin derivatives in the non-degenerate case we are dealing with. In this context, 
it is crucial to bound the first Malliavin derivative of Xf (called DXf in the sequel) effi- 
ciently This is where our asymmetry hypothesis on the vector fields Vi, . . . ,Vn enter into 
the picture, and we shall see (at Theorem 13. 2p how asymmetry properties yield an easy 
deterministic bound on DX^^. This result enables to get concentration results for the law 
of X^, and is the key to our density bounds as well. 

As another interesting consequences of the deterministic bound on DXf , we also obtain 
Log-Sobolev inequality and Poincare inequality for the law of Xt. 

Once the picture for the smooth case (when > |) becomes clear, we are able to 
extend some of our results described above to the irregular case when ^ < H < ^. In 
particular, we are able to prove 

Theorem 1.4. Fix G (|,i). A ssume Hypothesis Let X^ be the solution to 

equation ^ and '-/Xt ihe Malliavin matrix of X^ , t > 0. We have |det7xt|~"'^ G L°°(P). 
The random variable Xf admits a smooth density pxit, ■) and for any 6 < H there exist 
2 positive constants , 4 such that 

(4) px{t,y)<ci'^exp(^-ci'^\y\'y 
for ally G R"^. 

The existence of a density for solutions to stochastic differential equations of the form ([2]) 
under Hormander's condition has been obtained by Cass and Friz [6] for any | < < |. 
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While finishing the current article, an important step towards the study of regular densities 
in the rough case j < H < ^ has been accomplished in j8], where integrability estimates 
for the Jacobian of equation (|2]) are established. Nevertheless, as of today, besides the 
result of P. DriscoU [12], when H < ^, to the best of our knowledge. Hypothesis 11.21 is a 
first wide class of examples where we have an affirmative answer for the smoothness of 
the density. Our bound on the inverse of the Malliavin matrix together with polynomial 
bounds on the Holder norm of the Malliavin derivative allows then to obtain the sub- 
Gaussian upper bound (jl]). 

Notations: Throughout this paper, unless otherwise specified we use | ■ | for Euclidean 
norms and || ■ \\lp for the norm with respect to the underlying probability measure P. 

Consider a finite-dimensional vector space V. The space of valued Holder continuous 
functions defined on [0, 1], with Holder continuity exponent 7 G (0, 1), will be denoted by 
C^iy), or just when this does not yield any ambiguity. For a function g G C^iV) and 
0<s<t<l,we shall consider the semi-norms 

(5) llfl'IUt,7 = sup 

s<u<v<t 



\9v - 9u 


\v 




\v — u\ 


7 



The semi-norm ||5'||o,i,7 will simply be denoted by \\g\ 



7- 



2. Stochastic calculus for fractional Brownian motion 

For some fixed if G (|, 1), we consider (fi, J-", P) the canonical probability space associ- 
ated with the fractional Brownian motion (in short fBm) with Hurst parameter H. That 
is, Q = Co([0, 1]) is the Banach space of continuous functions vanishing at equipped 
with the supremum norm, J-" is the Borel sigma-algebra and P is the unique probability 
measure on Q such that the canonical process B = {Bf = [B], . . . , Bf), t G [0, 1]} is a 
fractional Brownian motion with Hurst parameter H. In this context, let us recall that 
-B is a (i-dimensional centered Gaussian process, whose covariance structure is induced by 
equation ([1]). This can be equivalently stated as 

R (t, s) := E [Bl B{] = ^ (s^^ + - |t - s\^^) , for s, t G [0, 1] and j = 1, . . . , d. 

In particular it can be shown, by a standard application of Kolmogorov's criterion, that 
B admits a continuous version whose paths are 7-Holder continuous for any 'j < H. 

This section is devoted to give the basic elements of stochastic calculus with respect to 
B which allow to understand the remainder of the paper. 

2.1. Malliavin calculus tools. Gaussian techniques are obviously essential in the anal- 
ysis of fBm, and we proceed here to introduce some of them (see [26] for further details): 
let S be the space of M'^- valued step functions on [0, 1], and "H the closure of S for the 
scalar product: 

d 

((l[0,ti], ■ ■ ■ , Mo,td])^ (l[0,^si], ■ ■ ■ , l[0,sd]))w = ^R{U, Si). 

1=1 
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We denote by K"^ the isometry between H and -^^^([0, 1]). When if > ^ it can be shown 
that 1],M'^) C n, and when | < iJ < ^ one has 

CHC L\[0,1]) 

for all 7 > I - if. 

Some isometry arguments allow to define the Wiener integral B{h) = {hs,dBs) for 
any element /i G "H, with the additional property E,[B (hi) B {h2)] = {hi, h2)n for any 
hi,h2 ^T-i- A J-"- measurable real valued random variable F is then said to be cylindrical 
if it can be written, for a given n > 1, as 

F = f {B{h'), B{h-)) = f{j\hl, dBs), j\h:, dBs)) , 

where E H and / : M" — )■ M is a C°° bounded function with bounded derivatives. The 
set of cylindrical random variables is denoted S. 

The Malliavin derivative is defined as follows: for F G 5, the derivative of F is the M'^ 
valued stochastic process (I)tF)o<t<i given by 

BtF = Y,h\t)^ {Bih'), . . . ,Bih^)) . 

i=l * 

More generally, we can introduce iterated derivatives. If F G iS, we set 

For any p > 1, it can be checked that the operator D'^ is closable from S into L^(f2; l-L®^). 
We denote by ©^'^("H) the closure of the class of cylindrical random variables with respect 
to the norm 

ll^ll.,P= (E(n + EE(||D^F||;,J 

V i=i 

and 

p>l k>l 

2.2. Differential equations driven by fBm. Recall that we consider the following kind 
of equation: 

(6) X^ = x+ f Vo{X^,)ds + Y.f V^i^s)dBl, 

Jo ^0 

where the vector fields Vq, . . . ,Vn are C°°-bounded. 

When equation Qj is driven by a fBm with Hurst parameter if > ^ it can be solved, 
thanks to a fixed point argument, with the stochastic integral interpreted in the (pathwise) 
Young sense (see e.g. [H]). Let us recall that Young's integral can be defined in the 
following way: 
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Proposition 2.1. Let f E , g E C with + k > 1, and < s < t < 1. Then 
the integral g^df^ is well-defined as limit of Riemann sums along partitions of [s,t]. 
Moreover, the following estimation is fulfilled: 



(7) 



9cdf( 



<C\\fUgUt-s\\ 



where the constant C only depends on 7 and n. A sharper estimate is also available: 



(8) 



<\9s\\\fh\t-s\^ + c,JfUg\Ut-s\^+^ 



With this definition in mind, we can solve our differential system of interest, and the 
following moments bounds are proven in |19) : 

Proposition 2.2. (Hu-Nualart) Consider equation ^ driven by a fBm B with Hurst 
parameter H > ^. Let us call its unique P-Holder continuous solution, for any (] < H. 
Then 

(1) When vector fields V are C°° -bounded, we have 

sup |Xf| < |x| + cv,t\\B\\]I^ p. 
te[o,T] 

(2) If we only assume that vector fields V have linear growth, with VV,V^l^, bounded, 
the following estimate holds true: 



(9) sup \X^\ < (1 + \x\) exp (cv,t\\B\ 



0,T,/3 



tG[0,T] 

Remark 2.3. The framework of fractional integrals is used in [H] in order to define integrals 
with respect to B. It is however easily seen to be equivalent to the Young setting we have 
chosen to work with. 

When the Hurst parameter | < if < |, equation ([6]) can be solved, again by fixed point 
argument, with the stochastic integral interpreted in the (pathwise) rough path theory 
(see e.g. jl^ and |22]). In this case, we obtain 



Proposition 2.4. (Besalu-Nualart Consider equation ^ driven by a fBm B with 
Hurst parameter ^ < H < ^. Denote by its unique ^-Holder continuous solution, for 
any P < H . If the vector fields V are -bounded, then for any A > and 6 < H 



E l^expA sup iXtl" j j < 00. 

Once equation (|6]) is solved, the vector X^ is a typical example of random variable 
which can be differentiated in the Malliavin sense. In fact, fix if G (|, 1), one gets the 
following results (see and |28] for further details): 

Proposition 2.5. Let X^ be the solution to equation ^ and suppose Vi 's are - bounded 
vector fields on . Then for every i = 1, . . . ,d, t > 0, and x G M"', we have X^'^ G D°°('H) 
and 

TDiX^ = Jo^tJo\sV,{X,), j = l,...,d, 0<s<t, 
where D^Xf'* is the j-th component o/D^Xf'*, and Jo~^t = 
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Finally the following approximation result, which can be found for instance in [13], will 
also be used in the sequel: 

Proposition 2.6. For m > 1 and T > 0, let B"^ = {-B™; t E [0,T]} be the sequence 
of linear interpolations of B along the dyadic subdivision of [0, T] of mesh m; that is if 
C = «2-™T fori = 0,...,2"*; then for t E 

51" = Bf"i -| (Btm — Bfrn) . 

t "-i j.m j-m \ "-i+i ''i ' 

Consider X™' the solution to equation ^ restricted to [0,T], where B has been replaced 
by S™. Set also D^.Xf = Jo^jJo^,V,(Xf ), for ] = I, . . . ,d and Q < s < t. Then almost 
surely, for any 7 < if and t E [0, T] the following holds true: 

(10) lim (||X^ - + llD^Xf - D^Xfl^) = 0. 



m— >-oo 



3. Estimates for solutions of SDEs driven by fBm: the smooth case 

Throughout this section, we fix H E (|, 1)- Recall that designates the solution 
to ([6]). This section is devoted to get some further bounds for Xf and its Malliavin 
derivatives, under Assumption ll.li 

Notice that among our set of hypothesis, the antisymmetric property ([3]) for the vector 
fields Vi, . . . ,Vn is the most specific one. It will be mainly used through the following 
lemma: 



Lemma 3.1. Let Ai,A2 be n x n matrices, whose exponential are defined by e^^ = 
Xl^o^i/^' /o?" J = 1,2. If we assume that is skew symmetric, then 

where \\A\\ stands for the Euclidean norm of a matrix A. 

Proof. Let us first prove the following (presumably classical) identity: 

(11) e*(^^+^^) = e*^^ - r e(*-^)(^^+^2Me^^^ds. 

Jo 

Indeed, consider the function s ^-)■ (f{s) defined on [0, 1] by Lp{s) = e(*~^)("^i+"^2)gsA2^ Then 
it is easily seen that (p is differentiable and 



By writing 







relation (fTT]) is now easily obtained. 

Let us see now the implications of ( ITT]) : according to the fact that A2 is skew-symmetric, 
we have ||e*^^|| < 1 for any s > 0. Therefore, 



t 



< 1 + p^ll f 

Jo 
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By denoting f{t) = lle*^^^"'"^^)!! -^g t^us get 

/(t)<l + Pi|| ['f{s)ds. 
Jo 

This implies f(t) < e'l^^H* by a standard application of Gronwall's lemma and finishes our 
proof. 

□ 

We are now ready to prove the main result of this section, which is an almost sure 
deterministic bound for the Malliavin derivative of X^: 

Theorem 3.2. Under Assumption the Malliavin derivative of the solution to 
equation ^ can be bounded as follows for any T e [0, 1] (in the almost sure sense): 

(12) IIDX^II^ < Mexp(CT), with M = sup sup 

xGK" ||A|j<l 

and where the constant C linearly depends on Vq. In particular, one also has ||DXy||^ < 
Mexp(Cr). 

Proof. Let us focus on the proof of (fT2|) . Indeed, since is dominated by the supremum 
norm when if > |, this will be sufficient in order to prove the second claim of our theorem. 
We now split our proof in two steps. 

Stepl: Matricial expression for the derivative. Let us first restate Proposition 12.51 in the 
following form: Y)X^ is solution to 

T>{X^ = Jo^t($:V,)(x), < s < T, 

where $*V^- denotes the pullback action of the diffeomorphism $s = : M'' — )■ on the 
vector field Vj. Now, a simple application of the change of variable formula for Young 
type integrals yields 

d 

d{^:v,){x) = {^:[Vo,v,]){x)ds +Y.^^:[v,,v,]){x)dBi 

i=l 

Moreover, recall that the Lie brackets [\^, Vj] can be decomposed, according to Assump- 
tion [TTTl into 

d d 

[^o,VS-] = 5^<-V^, and [y„y,] = ^4l-,, 

k=l k=l 

with oj'lj = —ojI^ for i,j,k> 1. Hence, 

d d d 

k=l k=l i=l 

By denoting Ms the d x d matrix with columns 

Mi = BiX^ = J,^ri^:V,){x), 



J2\^Viix) 



i=l 
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we therefore obtain the equation 

(13) dMs = Ms (uJo{X:)ds + Y,<X:)dB\ ) , Mt = V{X^), 



1=1 



where V{X^) is the matrix with columns Vj{X^), 1 < j < d and where uJi{X^) is the 
skew symmetric matrix with entries u^jlX^). 

Step2: Approximation procedure. In order to show that the process Ai is uniformly 
bounded, consider the dyadic approximation introduced at Proposition 12.61 By applying 
ffTOl) to the couple {X, A^), it is sufficient to prove our uniform bounds on D^XJp, uniformly 
in m. Let us thus consider A^™ the solution of (fT3|) where B is replaced by B"^, that is 

dMT = MT (^o{XT)ds + J2 M^DdBT'''^ , MT = V{X^). 
In the sequel set also 



k=i 



TDK JDK 

-Dfm iD-im 

ARi := — for 1< n < 2™ and 1< < rf. 

fm _ fm ^ — — — — 

Then, for s G [tT_i,tT), we have 

dMT = MT Lo{XT)ds + Y,MXT) ^B^-_,t^ds 

\ k=l 

Therefore, for s G we obtain 



MT = MTn.e 



n — 1 " / 



Proceeding inductively, we end up with the following identity, valid for t G [tT-i^tT) ^^'^ 
n = 0,...,2'": 



-( W0(X!")ds + y;f-, r^m™ Uli(XY')AB^m ,m ds] 

(14) MT = v{x:T)e v*2™-i ^ »^ 

Owing to the skew-symmetry of Uk ioi k > 1, we can now apply Lemma (HU] to expression 
(HM in order to get 

ll-^ni <exp (^||a;o(X,|_^)||) ••■exp((C-t)||u;o(X,'^_J||) II^^t)!! 

(15) <Mexp(Cr). 

This is our claimed uniform bound on M[", from which the end of our proof is easily 
deduced. 

□ 

Once the bound (fT2l) on ||DXf ||^ is obtained, one can also retrieve some information 
on the Holder norms of DXf improving the general estimate (IHl). This is the content of 
the following proposition: 
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Proposition 3.3. Consider | < 7 < and set Cq = Mexp(CT), which is the constant 
appearing in relation (Q^. Under Assumption M . 1\ the Malliavin derivative of the solution 
to equation ^ can be bounded as follows for any T G [0, 1]; 

(16) ilDX^II^ < CTy,d (1 + \x\ + 

for a strictly positive constant CTy,d- 

Proof. We have shown at Theorem 13.21 that DX^ is governed by equation (fT3l) , and that 
||A^||oo ^ Cq- We will now separate our proof into a local and a global estimate, and 
notice that the constants appearing in the computations below might change from line to 
line. 

Stepl: Local estimate. Consider < s < t < T and set e = t — s. Let n < be two 
generic elements of Applying relation ([8]) to the expression of M.^ — Aiu given by 

equation f fT3|) . we obtain 

\Mv — Mu\ < Cq Cv\v — u\ 

d 

+ {\Mu\\uJ^{x:)\\\BUv - + \\Muj,{x^)\Ut^4BUv - ^p^) . 

i=l 

We thus obtain, for a constant Cy depending on the vector fields V , 

\\M\\s,ta < clcv\t - s\^~^ + dclcv\\B\\^ + [coCi/||X^||t- + cv||A^|Ut,7] II^IM^ - 
(17) 

< clcv {T^~^ + d \\B\\y) +dcl^cv\\X''mB\\y\t - + dcv\\B\\^\t - s\'^\\M\\s,t,'y 

Take now t — s = e such that (icy||-B||-y£:'^ = 1/2, namely e = [2d cv\\B\\y]~^^"' . Recall also 

that ||X'^||7 < c(l + ||i?||y''^) according to ^3\. It is then easily seen that relation (fT7|) 
yields 

||-M|Ut,7 < CQCT,v,d (1 + \x\ + ll^liy'^) := aT,V,d,B, 
for a strictly positive constant CT,v,d- 

Step2: Global estimate. We consider now s,t G [0,T] such that ie < s < {i + l)e < je < 
t < {j + l)e, where e has been defined at Step 1. Set also tj = s, tk = ke for i + 1 < k < j, 
and tj^i = t. Then 



\Mt-Ms 



^ (^T,V,d,B — ikV 



(18) < aTy,dA3-^ + ^) \t-s 

where we have used the fact that r 1— )■ is a concave function. Note that the indices i, j 
above satisfy (j — i + 1) < 2T /e. Plugging this into the last series of inequalities, we end 
up with our claim ([TBI) . 

□ 
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3.1. Log-Sobolev inequality. In this section, we present some interesting functional in- 
equalities which are usually studied in a Markov setting; namely, the logarithmic-Sobolev 
inequality and Poincare inequality. As we will see, these inequalities become available in 
our non-Markov case when we have uniform boundedness for the Malliavin derivative of 
(see Theorem I3.2p . 

We start with the following version of logarithmic Sobolev inequality for the law of X^. 
Theorem 3.4. Let C and M be in Theorem\3M We have for all / G and T G [0, 1], 

E/(Xr)'ln/(XT)' - In (E/(Xt)2) < 2M^e^^^T^^¥.\S/ f{XT)\\ 

provided the right hand side in the above is finite. 

Proof. The proof is standard by applying Clark-Ocone formula (see e.g. j5]). First recall 
the representation of fractional Brownian motion 

ft 



Bt= [ KH{t,s)dWs. 
Jo 



Here is a d-dimensional Wiener process. Denote by the Malliavin derivative with 
respect to the Wiener process W. We have 

(19) K*hB = JD^, 

where K'^ is the isometry from "H, the reproducing kernel space of B, to L^. By Clark- 
Ocone formula we have 

/(Xt)-E/(Xt)= / E[BYf{XT)\^,]dWs= [ E[K*Hi-DfiXT))s\^s\dWs. 
Jo Jo 

Hence, if we denote M, = E(/(Xt)|^s), < s < T, we have 

dMs = E[KU'DfiXT))s\^s]dWs. 

For simplicity we may assume that / > e for some e > 0, which can be removed afterwards 
by letting e tend to 0. Applying Ito's formula to Mg In Mg, we get 

E/(Xt) ln(/(Xr)) - E/(Xt) In (E/(Xt)) = E(MTlnMT) - E(MolnMo) 

1„ 1 ,,.,2 



(20) = -Ej^ —\E[K*H{Bf{XT))s\^s]\ds. 

Replace now / by in the above. By the Cauchy-Schwarz inequality, 

\E[K*j,iBfiXTf)s\^s] r = 4 |E[/(Xr)(V/(XT), KU^Xt)s)\^s\ \' 



< 



4E [fiXTfl^g] E [{VfiXT),K*H{DXT)s)'\^s] ■ 



Substituting the above to ( l20l) . together with Theorem 13.21 we obtain the desired result. 

□ 

As a corollary of the logarithmic Sobolev inequality obtained above, we have the fol- 
lowing Poincare inequality (see e.g. |T8| Theorem 8.6.8]). 

Theorem 3.5. Let C and M be in Theorem\3M We have for all f e C^, 

EfiXrf - {Ef{XT)Y < M^e^^'^T^"E\VfiXT)\^ 
provided the right hand side in the above is finite. 
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Remark 3.6. In the above, assume further that the vector fields Vi,...,Vrf form an uni- 
form elhptic system, we obtain the following natural expression of logarithmic-Sobolev 
inequality and Poincare inequality when working on a Riemannian manifold 

d 

i=l 

and 

d 

1=1 

3.2. Concentration inequality. It is classical that the boundedness of the Malliavin 
derivative in "H implies the Gaussian concentration inequality, more precisely (see [3T| 
Theorem 9.1.1] or |21)): 

Lemma 3.7. Let F G 3^''^{'H) such that almost surely ||DF||^ < C, where C is a non 
random constant. Then, for every 6 > 0, 

and therefore for every x > 0, 

(21) P(F-E(F) > x) < e"^. 

As a corollary of this lemma, we deduce 

Proposition 3.8. Assume that the Assumption \1.1\ is satisfied, then there exist C and 
M such that for every T > and A > 0, 



sup \X^\ - E sup |Xf I ) > A ) < exp 

0<t<T Vo<t<T 







Proof. Let 



F = sup \X^ 

0<t<T 



X I 



By Theorem 13. 2[ it is not hard to see that (see e.g. |26| ) 

||DF||^ < Me^^T^, 

where M and C are the same as in Theorem 13.21 Now an easy application of Lemma 13.71 
completes our proof. □ 

Remark 3.9. Notice that this relation can only be obtained for if > |. Indeed, denote by 
C° the space of continuous functions endowed with the supremum norm. Then the norm 
involved in pT| Theorem 9.1.1] is ||DXj||j;^oo(^). This norm is dominated by ||DXt||j;^oo((jO) 
only if i7 > i. 
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3.3. Gaussian upper bound. One natural way to estimate tlie density of Xt is to apply 
the results in [26| Chapter 2]. More precisely, we first have the following integration by 
parts formula for non-degenerate random vectors. 

Proposition 3.10. Let F = be a non- degenerate random vector and ■-yp the 

Malliavin matrix of F. Let G G D°° and (p be a function in the space G^{M.'^). Then for 
any multi-index a G {1,2,..., d}'', k > 1, there exists an element Ha G such that 

E[dMF)G] = E[ip{F)Ha]. 
Moreover, the elements are recursively given by 

d 

Ha = Haf.{H(^ai,...,ak-i))y 

and for l<p<q<oowe have 

\\Ha\\LP < ^p,gll7F^-'-^-^llfc,2*-vll^lU>9' 

where - = - + -. 

p q r 

Remark 3.11. By the estimates for Ha above, one can conclude that there exist constants 
/3, 7 > 1 and integers m, n such that 



(22) 



l-f^alUp < C'p.gll det7p^||2'^(||DF||),I' IIGIIfc.q. 



Remark 3.12. In what follows, we use Ha{F,G) to emphasize its dependence on F and 
G. 

As a consequence of the above proposition, one has 

Proposition 3.13. Let F = {F^, F'^) be a non-degenerate random vector. Then the 
density pf{x) of F belongs to the Schwartz space, and 

Pi.(x) =E[I|^>,}i7(i,2,...,d)(F,l)]. 
Now we state and prove a global Gaussian upper bound for the density function of X, 



X 

t ■ 



Theorem 3.14. Denote by Px{t,y) the density function of . There exist positive 
constants c[^\ c[^^, C3, and C4 such that for all t G [0, 00), 



(23) 



Px{t,y) < cS^^exp 



I (2) 

y\ - cU 



Here cl^"^ is of the form 



.(1) 



0,1 



(2) 

for some positive number a; and q ^ converges to a constant as t J, 0. 



as t J, 0, 
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Proof. Fix P e !)• By fl22|) and Proposition I3.13[ we have 

(24) <C(P{Xi>y})i||det7x!lir^l|D^i||^,,, 

for some constants a, 7 > 1 and integers m, n. Without lose of generahty, we may assume 
Ui > for 1 < i < d. Since otherwise, for example yj < 0, we can consider the alternative 
expression for the density 

and deduce similar estimate. In the following, we provide estimate for ||DXt||^^, || det Txt^H^o 
and P{Xj > y} respectively. 

By Proposition 12.21 there exist constants C > 0, depending on V,x and k, such that 

(25) sup \Xl\ < \x'\ + CT\\BfJ^ 



0<t<T 



/3' 



(26) sup ||7x, II < 



C 



1//3 



0<t<T ■ T2^^'^ 
rk'" ri t 

On the other hand we have, for some constant Afg (cf. |19|): 



(27) sup \Bil...BilXi\<Ce^^^^^^^o,T,p, 

0<t,ri<T ^ 



(28) P{||fi||o,r,;3 > r} < M^e 2T^iH-p) _ 

Hence 

II det7^^^||™a < 00 and ||DXf||^^^ < 00. 

Moreover, by f l26l) , f l271) and the tail estimate (l28l) , there exists constant C > independent 
of t such that 

II det rxX^ < ^ (1 + t'^nt. - 1, C*)^) , 



and 



In the above 



|DX,||^^^ < C (1 + r^('=+^»+^F(t, - 1, 1/(5, C)^ 



F{t,a,b,M)= / r'^e 2*2(^^-/3) g*"^ (ir 

for any a, M > and 2 > 6 > 0. Now set 

C( — II aex 7jfj|^„ II jjyi-tll;,^^. 
Elementary computation shows, for some constant a > 

(29) d'^r^of^] ast;0. 



Finally we estimate F{Xt > y}. Define 

6= max{|X,i|,|X,2|,...,|Xf|}. 
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First note that 

since the Malhavin derivative |DX^*| < Me*"* for each i = 1, 2, ...d. Hence C,t has the same 
concentration property as specified in Proposition 13.81 

By (l25l) and the concentration property for (t, we conclude 



(30) P{.Y, > t,} < P.i6 > 1- < exp 



< 



exp 



V / 



V / 

Now f l24|) . fl29|) and fl30|) give us the desired upper bound for the density px{t, y)- □ 

Remark 3.15. There is also an upper bound for the constant Cj^'' in the above theorem as 
1 1 oo. Indeed, by some elementary computation one can show that 

c[^^ < t"e'^ as 1 1 oo, 

for some a, P > 0. 

4. Extension to the irregular case 

From now on, our purpose is to extend the previous to the case of a fBm with Hurst 
index | < if < |. This requires the introduction of some rough paths tools, which is 
the aim of the current section. We shall use in fact the language of algebraic integration 
theory, which is a variant of the rough paths theory introduced in [13] (we also refer to 
|15] for a detailed introduction of the topic). 



4.1. Increments. The extended pathwise integration we will deal with is based on the 
notion of increments, together with an elementary operator 5 acting on them. The al- 
gebraic structure they generate is described in [T^ [T5]. but here we present directly the 
definitions of interest for us, for sake of conciseness. First of all, for an arbitrary real 
number T > 0, a vector space V and an integer k > 1 we denote by Ck{V) the set of 
functions g : [0, T]'^ — )■ V such that Qti-tk ~ ^ whenever tj = tj+i for some i < k — 1. Such 
a function will be called a {k — l)-increment, and we set C^{V) = Uk>iCk{V). We can now 
define the announced elementary operator 6 on Ck{V): 

fc+i 

(31) S : C,{V) ^ C,^,{V), {Sg)t,..,,,, = Y.^-l)'-'9t....u-t,,,, 

i=l 

where ti means that this particular argument is omitted. A fundamental property of 6, 
which is easily verified, is that 66 = 0, where 66 is considered as an operator from CkiV) 
to Ck+2{V). We denote ZCk{V) = Ck{V) n Ker5 and BCk{V) = Ck{V) n lm6. 

Some simple examples of actions of 6, which will be the ones we will really use through- 
out the paper, are obtained by letting g E Ci and h G €2- Then, for any s,M,t G [0,T], 
we have 

(32) {^g)st = gt- gs, and {6h)sut = hst - - Kt- 

Furthermore, it is easily checked that ZCk{V) = BCk{V) for any k > 1. In particular, the 
following basic property holds: 
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Lemma 4.1. Let k >1 and h G ZCk+i{V). Then there exists a (non unique) f G Ck{V) 
such that h = 6f . 

Proof. This elementary proof is included in and will be omitted here. However, let 
us mention that fti...tk = (~l)''"^^^oti...ti, is a possible choice. 

□ 

Observe that Lemma l4.ll implies that all the elements h G C2(y) such that 6h = 
can be written a.s h = 6f for some (non unique) / G Ci{V). Thus we get a heuristic 
interpretation of S\c2{v)'- it measures how much a given 1-increment is far from being an 
exact increment of a function, i.e., a finite difference. 

Notice that our future discussions will mainly rely on fc-increments with k < 2, for 
which we will make some analytical assumptions. Namely, we measure the size of these 
increments by Holder norms defined in the following way: for / G C2{V) let 

(33) ||/L= sup and C^{V) = {f e C^iV); \\f\\, < oo} . 

Obviously, the usual Holder spaces C^iV) will be determined in the following way: for a 
continuous function g G Ci{V), we simply set 

(34) 11^11^= 115^11^, 

and we will say that g G Ci{V) iff \\g\\fj, is finite. Notice that || • ||^ is only a semi-norm on 
Ci{V), but we will generally work on spaces of the type 

(35) C^^,{V) = {g : [0,T] ^ V; go = a, \\g\\^ < oo} , 

for a given a & V, on which \\g\\^ defines a distance in the usual way. For h G C^^V) set 
in the same way 



(36) \\h\\-/,p = sup 



\h 



sut 



s,M,ie[0,T] 



\u - Spit - u\p 



inf <^ ||/ii||p.,/,_p.; h = Yhi, < p^ < p 



where the last infimum is taken over all sequences {hi G CslV)} such that h = hi and 
for all choices of the numbers pi G (0, fi). Then || ■ ||^ is easily seen to be a norm on Cs{V), 
and we set 

C^{V) := {heCsiV); \\h\\^<oo}. 

Eventually, let Cl~^{V) = Ufj^^iCt^iV), and notice that the same kind of norms can be 
considered on the spaces ZC^lV), leading to the definition of some spaces ZC^{y) and 

With these notations in mind the following proposition is a basic result, which belongs 
to the core of our approach to pathwise integration. Its proof may be found in a simple 
form in |15] . 
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Proposition 4.2 (The A-map). There exists a unique linear map A : ZC^'^iV) — t- C2^{y) 
such that 

In other words, for any h G C\^{y) such that 6h = there exists a unique g = A(/i) G 
Cl'^{V) such that 6g = h. Furthermore, for any fi > 1, the map A is continuous from 
ZC^iV) to Cl^iV) and we have 

(37) ||A/.||,<^^||/.|U, heZC^,{V). 

Let us mention at this point a first link between the structures we have introduced so 
far and the problem of integration of irregular functions. 

Corollary 4.3. For any 1-increment g E C2{V) such that 5g E Cl'^ , set Sf = {Id— AS)g. 
Then 

n-l 
i=0 

where the limit is over any partition Ust = {to = s, . . . ,tn = t} of [s, t], whose mesh tends 
to zero. Thus, the 1-increment Sf is the indefinite integral of the 1-increment g. 

4.2. Computations in C*. Let us specialize now to the case V = M, and just write 
for C2^(M). Then (C*, 6) can be endowed with the following product: for g E Cn and h E Cm 
let gh be the element of Cn+m-i defined by 

(38) (5'^)tl,...,W„+l = fi'tl,...,in^tn,-,im+n-l' ^ 1 , • • • , Wn- 1 G [0 , T] . 

In this context, we have the following useful properties. 

Proposition 4.4. The following differentiation rules hold true: 

(1) Let g E Ci and h E Ci. Then gh E C\ and 6{gh) = 6g h + g Sh. 

(2) Let g E Ci and h E €2- Then gh E C2 and 6{gh) = —6g h + g 6h. 

(3) Let g E C2 and h E Ci. Then gh E C2 and 6{gh) = 6g h + g dh. 

The iterated integrals of smooth functions on [0,T] are obviously particular cases of 
elements of C, which will be of interest for us. Let us recall some basic rules for these 
objects: consider / E Cf^ and g E C^, where Cf^ denotes the set of smooth functions on 
[0,T]. Then the integral / f dg, which will be denoted indistinctly by J f dg or J{f dg), 
can be considered as an element of Namely, let S2,t denote the simplex {(s,t) E 
[0, T]2 : < s < t < T}, for (s, t) E S2,t we set 



JMdg) = 




The multiple integrals can also be defined in the following way: given a smooth element 
h E and {s,t) E S2,t, we set 



Jst{hdg) = 
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In particular, for e , p G and p e the double integral Jstifdfdf) is 
defined as 

Jstif dfdf) =(^Jf dfdf^ = Jsu if dp) dfl 

Now suppose that the nth order iterated integral of f"'^^df"' ■ ■ ■ df"^, which is denoted by 
■ ■ ■ dp), has been defined for p G C^. Then, if p e C^, we set 

(39) JstiP^'dP ■ ■ ■ dpdp) = {r^'dP ■ ■ ■ dp) dp, 

which recursively defines the iterated integrals of smooth functions. Observe that an nth 
order integral J{df'^ ■ ■ ■ dpdf^) can be defined along the same lines, starting with 

J{dp = 6f, Mdp dp) = f Jsuidp) dp = f {5p)^^ dp, 

J s J s 

and so on. 

The following relations between multiple integrals and the operator 5 will also be useful. 
The reader is sent to [15] for its elementary proof. 

Proposition 4.5. Let f E and g E . Then it holds that 

5g = J{dg), 5 {J{fdg)) = 0, 5 {J{dfdg)) = {5p{6g) = J{dpJ{dg), 

and 

n-l 

s {j{dr ■ ■ ■ dp)) = Y,j{dr--- dp+') J {dp ■ ■ ■ dp) . 

i=l 

4.3. Weakly controlled processes. The rough path theory allows to define and solve 
differential equations driven by a generic Holder continuous path B provided enough 
iterated integrals of this function can be defined. We shall briefly recall how this is done, 
in the simplest nontrivial case of a Holder continuity exponent | < 7 < |- Observe that 
we keep here the notation B for the underlying path as in the fBm case for notational 
sake, while the theory can be applied to much more general situations. 

The basic assumption one has to add in order to define our objects when 7 > | is the 
existence of an (abstract) double iterated integral of B with respect to itself, which can 
be defined as follows: 

Hypothesis 4.6. The path B is W^-valued --y-Holder with 7 > | and admits a Levy area, 
that is a process e Cl^iW^''^) satisfying 

= ® B\ I. e. ^ = [B^'%^[B^'^U 

for s,u,t E iSs^T dn-d i,j E {!,..., rf}. We also assume that can be obtained in the 
following way: consider the sequence of linear dyadic approximation 5™ of B defined 
like in Proposition \2.(A For < s < t < 1 and 11,12 E {!,..., c/}, set _ 
X!<«i<«2<t ^-^u'P dB^^^"^ , which is defined as a Riemann-Stieljes integral. Then we suppose 
that B^'™ converges to B^ in the norm of . 

It should be noticed at this point that fBm satisfies the above assumption, as shown 
in P [131122]: 
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Proposition 4.7. Let B be a d- dimensional fBm with H > ^ as defined in Section 
For the linear dyadic approximation 5™ of B defined at Proposition \2.6[ the increment 
B^'"* almost surely converges to an element satisfying Hypothesis \4 ■ 6] The convergence 
holds in any norm for ^ < H . 

The first difference between tlie Young case and tlie situation of a Holder continuity 
exponent | < 7 < | is tliat a restriction lias to be imposed on the class of allowed 
integrands with respect to B. This class is called the class of weakly controlled processes, 
and is defined as follows: 

Definition 4.8. Let z he a process in C7(M") with | < 7 < | (that is, N := [I/7J = 2j. 
We say that z is a weakly controlled path based on B and starting from a if zq = a, which 
is a given initial condition in M", and Sz G (ffi") can be decomposed into 

(40) 6z' = rB^'''+r\ I.e. (5^%^ = eB^f + 

for all {s,t) G S2^t. In the previous formula, we assume ( G C7(ffi"''^), and r is a regular 
part such that r G C2'^{M.'^). The space of weakly controlled paths starting from a will be 
denoted by Q^^a(R"'), and a process z G Q^^a(IR") can be considered in fact as a couple 
{z,(). The natural semi-norm on Q^^a(]R") is given by 

^^[z; Q^,,(M")] = Ar[z;C7(M")] + M[C;C'{°{W'''')] +M[C;C^{W''')] + N'[r;C^\W)], 

withJ\f[g;C^] defined by ^ and Af[C;C^{V)] = supo<,<T |Cs|y- 

Two basic steps in order to define and solve differential equations with respect to B are 
then: 

(1) Study the decomposition of f{z) as weakly controlled process, when / is a smooth 
function and z a weakly controlled process. 

(2) Define rigorously the integral J z^dB^ = J{zdB) for a weakly controlled path z 
and compute its decomposition ( 140|) . 

We shall now detail a little this program. 

Let us see then how to decompose f{z) as a controlled process when / is a smooth 
enough function, a step for which we first introduce a convention which will hold true until 
the end of the paper: for any smooth function / : M" — )■ M, > 1, (ii, . . . , 2^) G {1, . . . , d}^ 

and ^ G M" , we set 

With this notation in hand, our decomposition result is the following: 

Proposition 4.9. Let f : W ^ be a Cl function such that f{a) = a, z a controlled 
process as in Definition \4.8\ and set z = f{z). Then z G Q-y^a{^), o-nd it can be decomposed 
into Sz = C'lB-*-'*^ + f, with 

C = dj{z) and f = [5f{z) - dj{z) 5z^] + dJiz) r\ 

Furthermore, 

(42) U[z- Q,,m] < Cf,T (1 + ^f'[z; QjA^n]) ■ 
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Let us now turn to the integration of weakly controlled paths, which is summarized in 
the following theorem. 



Theorem 4.10. For a given | < 7 < |, let B be a process satisfying Hypothesis \4-6 
Furthermore, let m G Q^^blMf^) with rriQ = 6 G M'' and decomposition 

(43) 6m' = fi''' B^''' + r\ where fi G CJiR"^'"^), r G C^'^iR'^). 
Define z by Zq = a ^R and 

(44) 6z = m' B^ ' + fi''^ B^ '^' - A (r* B^'* + 5^"^ B^''^') . 
Finally, set 

JstimdB) = / (m„, dBu)^d = 6zst- 

J s 

Then: 

(1) z is well-defined as an element of Q^^a{^), o-nd coincides with the Riemann- Stieltjes 
integral of m with respect to B whenever these two functions are smooth. 

(2) The semi-norm of z in Q^^a(]R) can be estimated as 

(45) Af[z- Q^,a(R)] <Cb{1+ Af[m; Q^,b(R'')]) , 
for a positive constant Cb which can be bounded as follows: 

Cb<c (U[B'^] C1{R'^)] + A/'[B^; C2^(M'^')] j , for a universal constant c. 

(3) It holds 

n-l 

(46) Mm dx) = hm ^ Yl K ^^'-^ + ^ ^5"- 



for any < s < t < T, where the limit is taken over all partitions Ugt = {s = to, . . . ,tn 
t} of [s,t], as the mesh of the partition goes to zero. 



4.4. Rough differential equations. Recall that we are concerned with equations of the 
form ([2]). In our algebraic setting, we will rephrase this as follows: we will say that is 
a solution to ([2]) if = x, G Q'y,x{R'^) and for any 0<s<t<lwe have 

(47) 6X:, = Vo{X:)du + J,, {V,{X^) dB') , 

where the integral J{yi{X^) dB') has to be understood in the sense of Proposition I4.1UI 
The following existence and uniqueness result is then classical in rough paths theory. 

Theorem 4.11. Let B be a process satisfying Hypothesis \4.6\ and Vq, ■ . ■ ,Vn a collection 
of vector fields which fulfill Hypothesis Then 

(i) Equation admits a unique solution X^ in Q'^^xiR'^)- 

(ii) Consider the linear approximation B"^ of B introduced in Hypothesis \4.6[ and set 
for the solution of the (ordinary) differential equation driven by the piecewise smooth 
function i?'". Then X"^ converges to X^ in norm. 
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(iii) Consider the sequence of dyadic partitions of Proposition \2. (A Define a process 
on the points t^ by Xq = x and 

V ( X^ ) 

(48) dX^ , = °^ + Vi (X" ) B,''' + Vi, V^, (X" ) B.^'f ^ . 

Complete the definition of X" on [0, 1] by linear interpolation. Then as n ^ oo, the 
process X" converges to X^ in Cj norm. 



Proof. We refer to [22] for the proof of the existence and uniqueness part, as well as to [T] 
for the same result in the algebraic integration setting. Part (ii) of our proposition stems 
from the continuity of the Ito map, which is also stated and proved in both fF2\ [Ti] . The 
approximation statement (iii) has first been stated by Davie [TU] and then been generalized 
in pS]. 

In the sequel we will simply try to relate the decomposition of the solution to Equa- 
tion (147|) as a controlled process and the numerical scheme given by (HHl) . a relation which 
turns out to be useful in the sequel. For this, we shall denote by r any increment in 
and by r" any increment in in the computations below, independently of their values. 
Observe then that, according to the right hand side of (147|) . the decomposition of X^' as 
a controlled process is given by 

SX""'^ = C^'JiB^'^'^ + r, with = Vl (X^). 

Hence, owing to Proposition 14. 9[ one has 

(49) SVi{X^) = C^^B^'^' + r\ with = djViiX'')^^' = V^^ViiX''). 

Now, if one desires an expansion of 6X^ up to increments of regularity l"*", consider 
again the right hand side of f H7|) . and compute it by a direct application of Theorem 14.101 
This yields 

6X:, = Vo{X:){t-s) + V,{X:)Bll' + C'^B^r + r^ 

1'* _L T/. T/r ^^■^■R^jii tt 



vo{x:){t -s) + mx:)Bti' + v,mx:)b 



St 



r" 



Thanks to identity (1461) . it is now easily seen that (HHj) is a natural candidate for our 
numerical scheme. 

□ 

We show now how to get efficient bounds on the solution to equation ( 147|) out of its 
numerical scheme. This step is understood as a warmup for the same kind of estimates 
concerning the Malliavin derivative of the solution. 



Proposition 4.12. Under the assumptions of Theorem \4.11\ the solution to equation 
satisfies: 

(50) 11X^1^ <cy(l + \\B%/^ + IIB^II^/'^' 

where Cy is a constant which only depends on the vector fields Vq, . . . , V^. 

Remark 4.13. This proposition is shown in by identifying the signature of B with 
the signature of a certain finite variation process, plus some easy estimates for ordinary 
differential equations. We have included here a direct elementary proof of (130]) because 
we haven't been able to find them in the literature under this form, and mostly because 
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the generalization of our estimates to linear cases will be obvious from the considerations 
below. 

Proof of Proposition \4-i^ Theorem 14. 1 1 1 part (iii) asserts the convergence of the approx- 
imation towards as — )■ oo. It is thus sufficient to prove relation (l50l) for X", 
uniformly in n. One can also be easily reduced to prove 

I A I 

-P ^;r^<-v i + l|B^ie + l|B^lir 

0<i<j<2" |tj ~ I ' ^ 

which is what we shall proceed to do. We now divide our proof in several steps. 
Step 1: Expression for dX'"" . Set 

(51) ql = Vo{X:) it-s) + V,{X:)Bll' + V,,V,,{X:)B',r\ 

so that SX^n^n = Xl«=i For i < j, we also construct a dyadic partition {rf; < 

k < K,0 < I < 2*''} of the set {t^, . . . ,t^} inductively in the following way: set Tq = ti 
and = tj. Now, if the rj'^s are known, we set T2i^^ = rj^- Furthermore, if r^i^^ = tm 
and T2i^\ = t^i, then take Tg;^^ = t^*, with m* = [(m + m')/2\. This procedure is then 
non trivial as long as j — i > 2'^, which means that we stop at K = [log2(j — i)] ■ Here is 
then a simple example of construction: consider i = 1, j = A. Then we have: 

r° = 1, r? = 4; = 1, = 2, = 4; = 1, = 1, = 2, = 3, = 4. 



With these notations in hand, it is easily checked that the relation 5^^" = Y2''i=i Qt^t" 



K-l K-l K-l. 

21 '21+1 '21+2 

1=0 1=0 



K 2'=-! 



» 3 

can also be written as 6X^n^n = J2k=o^ C^r^ • Furthermore, 

q!!^k^k + Q'l^fc = q^k-i k-l — (.Sq"')^k k k , 

^21^21+1 ^21+1^21+2 n + 1 ^ ^ '^21^21+1^21+2' 

and summing this equality for k = K and / = 0, . . . , 2^"^ — 1 we get 

2^-1 2^-1 2^-1 

1= 

Iterating, we obtain 

^ ' ^3 S' ^ ^ ^21^21 + 1^21+2 

k=l 1=0 

Step 2: Expression for 5q"'. Denote by / the identity function on R, so that Slgt = t — s. 
Start then from expression (l52l) and use Proposition 14.41 in order to get, for any s,u,t in 
the dyadic partition, 

dq:,, = 6Vo{X-)su {t-u) + 6V,{X-),^Bl:i + 6 [V,,V,,{X-)l^Blr ' V,,V,,{X:)6B',Cr , 
or otherwise stated thanks to convention (l38l) . 

(53) <5g" = 6Vo{X") 61 + 6V,{X'')B^'' + 6 [Vi^Vi^{X'')] B^''^'^ - V,yi^{X''')6B^''''\ 
Observe now that one can prove, as for ( 149|) . that 

5\/.(X") = Vi,V,{X'')B^''' + r''"', with |r^f I < cv\SX^t\^. 
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Thus, according to the fact that (5B^'*^*^ = B-'-'^^B-'-'*^, one can recast f l53|) into 

(54) = SVoiX") SI + r''"B^'* + S [V^,V^,iX'')] B^-'^*^ := p^'" + p^-- + p3,n_ 

Notice that for j = 1,2,3, the increment p^'" hes into Cp . Furthermore, it is readily 
checked that 

(55) Ip^;",! < cv\SX:J \t-u\'+\ IplZl < cv\6X:J'\\B^\\,\t-u\\ 
and \pl'Z\<cv\5X:J\\B%,\t~u\'^. 

Step 3: An induction procedure. Let us consider an integer £ > 1 and the quantity: 

I A I 
K= sup 

We locahze now our study to an interval of the form [a, a + 77] with an arbitrary positive 
number a, and 77 small enough. We will prove that if 77 is of order (1 + ||Bi||y^^ + 
l|B^||2(^^)-\ then < cy(l + ||Bi||^ + ||B2||2/^) by induction. 

The case i = 1 being trivial, let us assume that the hypothesis is true up to a given 
i > 1. Take now 1 < i < i and j = £ + 1. According to ( !52|) . write 

K 2'---l 

k=l 1=0 

In the right hand side of this decomposition, all the points T|;,r2;_,_i are of the form tp 
with p < i. Thus fl5^ . our bounds fl55|) on p^'" and the induction hypothesis entail 



I j.n. 1 7 

IS' I 

Hence, putting together the last two inequalities, taking into account that we work on 
an interval of size 77 and that we have chosen j = ^ + 1, we end up with the following 
induction relation: Ni^i < Fri{Ni), where the function : ]R_|_ — t- is defined by 

F.iO = e V [(1 + IIB^II^ + ||B2||2^ r]^) + (l + HB^H^^^ + l|B^||^e') v'^] ■ 

By separating the cases 77"^^ < 1 and rj'^^ > 1, one can also prove that -Fjy(0 — ^ ^ V9^(^), 
with 

<^,(0 = Cv {\\B%,r]'^ + \\B%,r]^^) + (l + HB^H^ + HB^Ib^r/^) := a^' + c. 

In order to obtain a bound of the form Ne < M which remains valid for all i, it is now 
sufficient to have the interval [0, M] left invariant by ifrj- 

We let the reader check the following elementary fact: whenever 4ac is of order 1, the 
interval [0, M] is left invariant by the application 1— )■ aC.'^+c, with M of order c. Applying 
this to ifj^, we obtain that 

(56) V ' - 



Furthermore, it is obvious from (ISTl) that 

<cv {l + \\B% + \\B%,\t]-t^\^) 



l + ||Bi||, + ||B2||^/^ 



iV,< llB^II.+ llB^ll^f 
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Step 4: Conclusion. We have thus obtained that on any interval of length rj given by fl56|) . 
we have < IIB"""!!^ + ||B^||2^|,^. Our claim (150|) is now easily deduced by dividing an 

arbitrary interval [s,t] into subintervals of length rj as in ( ITSl) . 

□ 

4.5. Estimates for the Malliavin derivative. We are now interested in extending 
Proposition 13.31 beyond the Young setting. Recall thus that we are concerned with equa- 
tion (fT3|) . which can be written in our algebraic integration setting as 

(57) 6Mst= I MuUJoiX:)du + Jst{MuiiX^)dB'), 0<s<t<T, 

J s 

with final condition Air = V{^t)- Then we have the following equivalent of Theorem 

EH 



Proposition 4.14. Theorem \4.11\ holds true for equation [5l\) under Hypothesis li.il 



and \4.6\ The discretization scheme for Ai can be written as: 
{51 



tfcti--4-1 tl, 



2^ 



Proof. As in the case of X^, we only justify expression (l58l) . Note that, according to 
equation (1571) . we have 6Aist = Cl^lt with (I = Aig^^iiX^). Hence, invoking Propo- 
sition HiH and relation ( l49l) . we get 

6[M uj,,] = M h,a;.,(X^) + V,,uj,,{X^)] B^ '^ + r. 

Expanding now the right hand side of equation (1571) with the help of Theorem 14.101 one 
easily gets 



5M=M [uJo{X'-) 61 + u,{X^)B^'' + {u,,u,,{X'^) + V,,u,,{X^)) B^-^^^^j + r" 



4 

which gives the desired justification of our scheme (l58l) . 

□ 

We are now ready to state and prove our bounds for the process M.: 

Proposition 4.15. Let M. he the unique solution to ( f57|) under Hypothesis and \4.6[ 
Then 

(i) The bound flE) on ||A^||oo still holds true in our irregular context. 

(ii) ||A^||-Y satisfies inequality (E^. 

Proof. By the continuity property of Ito's map for Ai, the approximation procedure of 
Ai by Ai"^ described in the proof of Theorem 13.21 is still valid. The desired bound for 
||A^||oo is thus obtained just like in (|T5|) . 

Once a bound on ||A^||oo is available, ||A^||7 can be bounded by considering the Davie 
type scheme (l58i) . along the same lines as for X^. Details are left to the reader for sake 
of conciseness. 

□ 
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4.6. Density upper bound. We finish this section by extending the density estimate 
and functional inequalities obtained in the smooth case (when > |) to the irregular case 
(when I < < |). We first show that in the rough case, we can obtain the smoothness 
of density of Xt under Hypothesis 11.21 Indeed, for this purpose, we only need to show the 
following integrability of Malliavin matrix. 

Theorem 4.16. Fix H A ssume Hypothesis \l.iA Let 7^ he the Malliavin matrix 

of Xf, we have \ det7x|~"'^ G L°°(P). Therefore Xf admits a smooth density. 

Proof. It suffices to show that there exists C > such that for all f G M"^, 
Recall that 

and that the d x d matrix Ai is such that its j-th column M.^ is given by 

Ml = BIX,. 

We have 

d 

>CV / {v^M{fds = C v^MsMjvds>C mm Af^iM^fvds. 
~tJo Jo Jo 

In the above, we used that "H C L'^{[0, 1]) for the second step, and A^"* is that described 
in Theorem 13.21 

On the other hand one can show, according to the description of A^™ in Theorem 13.21 

uniformly for some constant C > 0, if we assume Hypothesis 11.21 The proof is therefore 
completed. □ 

Another consequence of the above boundedness of the Malliavin Matrix 7^ is the fol- 
lowing. 

Theorem 4.17. Fix if G (|, |). A ssume Hypothesis Let px{t,y) denote the density 
function of Xf . There exist constants c^^'^ and cf^ such that 

Px{t,y) < cf^exp (-cf)|y|^) y G 

for any 6 < H . 

Proof. The proof is similar to that of Theorem 13.141 We apply Theorem 14.151 and Theo- 
rem l4.16l to obtain bounds for ||DX(||^^ and || det 'JxlW^a. The tail estimate for F{Xt > y} 
is derived by Proposition 12. 4[ The rest of the proof is then clear. □ 
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Remark 4.18. We are not able to obtain log-Sobolev inequality as in Theorem 13.41 in the 
rough case, since when H < ^ the Hilbert norm "H is not controlled by L°°([0,1]). On 
the other hand, by reproducing the proof in Theorem 13.41 together with the following 
interpolation inequality 

||-|k<C(||-||oo+||-||7), 

and estimates in Proposition 14.151 we are able to prove the following version of a log- 
Sobolev type inequality: for T G [0, 1] we have 

EfiXTflnfiXTf - [EfiXTf] In (E/(Xt)^) < C,,t (E|V/(X^)r^)'/^ 
For all p > 1. Here Cp^T is a universal constant independent of /. 
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